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Abstract. 

We present a comprehensive hybrid library of synthetic 
stellar spectra based on three original grids of model atmo¬ 
sphere spectra by Kurucz (1995), Finks et al. (1994), and 
Bessell et al. (1989, 1991), respectively. The combined library 
has been intended for multiple-purpose synthetic photometry 
applications and was constructed according to the precepts 
adopted by Buser & Kurucz (1992): (i) to cover the largest 
possible ranges in stellar parameters (Ten, log g, and [M/H]); 
(ii) to provide flux spectra with useful resolution on the uni¬ 
form grid of wavelengths adopted by Kurucz (1995); and (hi) to 
provide synthetic broad-band colors which are highly realistic 
for the largest possible parameter and wavelength ranges. 

Because the most astrophysically relevant step consists in 
establishing a realistic library, the corresponding color calibra¬ 
tion is described in some detail. Basically, for each value of 
the effective temperature and for each wavelength, we calcu¬ 
late the correction function that must be applied to a (theo¬ 
retical) solar-abundance model flux spectrum in order for this 
to yield synthetic UBVRIJHKL colors matching the (empiri¬ 
cal) color-temperature calibrations derived from observations. 
In this way, the most important systematic differences exist¬ 
ing between the original model spectra and the observations 
can indeed be eliminated. On the other hand, synthetic UBV 
and Washington ultraviolet excesses 5(u-b) and S(c-m) and 
5(cr-Ti) obtained from the original giant and dwarf model spec¬ 
tra are in excellent accord with empirical metal-abundance cal¬ 
ibrations (Lejeune & Buser 1996). Therefore, the calibration 
algorithm is designed in such a way as to preserve the orig¬ 
inal differential grid properties implied by metallicity and/or 
luminosity changes in the new library, if the above correction 
function for a solar-abundance model of a given effective tem¬ 
perature is also applied to models of the same temperature but 
different chemical compositions [M/H] and/or surface gravities 
log g. 

While the new library constitutes a first-order approxima¬ 
tion to the program set out above, it will be allowed to develop 
toward the more ambitious goal of matching the full require¬ 
ments imposed on a standard library. Major input for refine¬ 
ment and completion is expected from the extensive tests now 
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being made in population and evolutionary synthesis studies 
of the integrated light of globular clusters (Lejeune 1997) and 
galaxies (Bruzual et al. 1997). 


1. Introduction 

The success of population and evolutionary synthesis calcula¬ 
tions of the integrated light of clusters and galaxies critically 
depends on the availability of a suitable library of stellar spec¬ 
tral energy distributions (SEDs), which we shall henceforth call 
stellar library. Because of the complex nature of the subject, 
there are very many ways in which such calculations can con¬ 
tribute to the solution of any particular question relevant to the 
stellar populations and their evolution in clusters and galaxies. 
In previous studies, the particularities of these questions have 
largely determined the properties that the corresponding stel¬ 
lar library must have in order to be considered suitable for the 
purpose. 

There is now a considerable arsenal of observed stellar li¬ 
braries (for a recent compilation, see e.g. Leitherer et al. 1996), 
each of which has its particular resolution, coverage, and range 
of wavelengths as well as its particular coverage and range of 
stellar parameters - but which, even if taken in the aggre¬ 
gate, fall short by far of providing the uniform, homogeneous, 
and complete stellar library which is required now for a more 
systematic and penetrating exploitation of photometric popu¬ 
lation and evolutionary synthesis. 

Ultimately for this purpose, what is needed is a uniform, 
homogeneous, and complete theoretical stellar library, provid¬ 
ing SEDs in terms of physical parameters consistent with em¬ 
pirical calibrations at all accessible wavelengths. Thus, the 
above goal can be approached by merging existing grids of 
theoretical model-atmosphere spectra into the desired uniform 
and complete stellar library, and making it both homogeneous 
and realistic by empirical calibration. 

A variant of this approach was first tried by Buser & Ku¬ 
rucz (1992), who constructed a more complete theoretical stel¬ 
lar library for O through K stars by merging the O-G-star 
grids of Kurucz (1979a,b) with the grids of Gustafsson et al. 
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(1975), Bell et al. (1976), and Eriksson et al. (1979) for F- 
K stars. In their paper, Buser & Kurucz solved for unifor¬ 
mity and homogeneity by recomputing new late-type spectra 
for Kurucz’s (1979) standard grid of wavelengths and using 
the Kurucz & Peytremann (1975l-|atomic opacity source ta¬ 
bles. The resulting hybrid libraryu has, indeed, significantly 
expanded the ranges of stellar parameters and wavelengths for 
which synthetic photometry can be obtained with useful sys¬ 
tematic accuracy and consistent with essential empirical effec¬ 
tive temperature and metallicity calibrations (Buser & Fenkart 
1990, Buser & Kurucz 1992, Lejeune & Buser 1996). 

In the meantime, Kurucz (1992, 1995) has provided a 
highly comprehensive library of theoretical stellar SEDs which 
is homogeneously based on the single extended grid of model 
atmospheres for O to late-K stars calculated from the latest 
version of his ATLAS code and using his recent multi-million 
atomic and molecular line lists. The new Kurucz grid - as we 
shall call it henceforth - indeed goes a long way toward the 
complete library matching the basic requirements imposed by 
synthetic photometry studies in population and evolutionary 
synthesis. As summarized in Table 1 below, SEDs are provided 
for uniform grids of wavelengths and stellar parameters with 
almost complete coverage of their observed ranges ! These data 
have already been widely used by the astronomical community, 
and they will doubtlessly continue to prove an indispensable 
database for population and evolutionary synthesis work for 
years to come. 

In the present work, we shall endeavor to provide yet an¬ 
other indispensable step toward a more complete stellar library 
by extending the new Kurucz grid to cooler temperatures. This 
extension is particularly important for the synthesis of old stel¬ 
lar populations, where cool giants and supergiants may con¬ 
tribute a considerable fraction of the total integrated light. 
Because model atmospheres and flux spectra for such stars - 
the M stars - have been specifically calculated by Bessell et al. 
(1989, 1991) and by Fluks et al. (1994), our task will mainly 
be to combine these with the new Kurucz grid by transforma¬ 
tion to the same uniform set of wavelengths, and to submit the 
resulting library to extensive tests for its realism. In fact, as 
shall be shown below, the process will provide a complete grid 
of SEDs which is homogeneously and consistently calibrated 
against observed colors at most accessible wavelengths. 

In Sect. 2, we shall briefly describe the different libraries 
used in this paper and the main problems that they pose to 
their unification. Because the spectra exhibit systematic dif¬ 
ferences both between their parent libraries and relative to 
observations, we set up, in Sect. 3, the basic empirical color- 
temperature relations to be used for uniformly calibrating the 
library spectra in a wide range of broad-band colors. This cal¬ 
ibration process is driven by a computer algorithm developed 
and described in Sect. 4. The actual color-temperature rela¬ 
tions obtained from the corrected library spectra are discussed 
in Sect. 5, and the final organization of the library grid is pre¬ 
sented in Sect. 6. In the concluding Sect. 7, we summarize 
the present state of this work toward the intended standard li- 


^The Buser-Kurucz library is a hybrid library in the sense that 
it is based on two distinct grids of model atmospheres which 
were calculated using different codes (ATLAS AND MARCS, 
respectively); it is quasi-homogeneous, however, because for 
both grids of model atmospheres the spectrum calculations 
were obtained using the same opacity source tables. 


brary, and we briefly mention those necessary steps which are 
currently in process to this end. 

2. The basic stellar libraries 

The different libraries used are from Kurucz (1995), Bessell 
et al. (1989, 1991), and Fluks et al. (1994) - which we shall 
henceforth call the K-, B-, and F-libraries, respectively. Al¬ 
though the K-library covers a very wide temperature range 
(from Tefi=50,000 K to 3500 K), it does not extend to the very 
low temperatures required to model cool AGB stars. These 
stars are very important for population and evolutionary syn¬ 
thesis, since they can represent up to 40% of the bolometric 
and even up to 60% of the K-band luminosity of a single stel¬ 
lar generation (Bruzual & Chariot 1993). It is natural, then, to 
provide the necessary supplement by employing the suitable li¬ 
braries that were specifically calculated for M-giant stars in the 
temperature range 3800-2500 K by Bessell et al. (1989,1991) 
and by Fluks et al. (1994). 

Table 1 summarizes the coverage of parameters and wave¬ 
lengths provided by these three libraries, and Figs. 1 and 2 
illustrate original sample spectra as functions of metallicity for 
two temperatures. 



Fig. 1. K—library dwarf model spectra [logg = 5) for a range in 
metallicity. All the spectra are normalized at A = 817nm. 

We should note the following: first, while both the K- and 
the B-library spectra are given for different but overlapping 
ranges in metallicity, the F-library spectra have been available 
for solar abundances only; secondly, the B-spectra are given for 
wavelengths AA > 491nm and nonuniform sampling, while the 
F-spectra also cover the ultraviolet wavelengths AA > 99nm 
in a uniform manner. 

As it is one of the goals of the present effort to eventu¬ 
ally also allow the synthesis of the metallicity-sensitive near¬ 
ultraviolet colors (e.g., the Johnson U and B), the B-library 
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Table 1. Parameter and wavelength coverage of the different libraries 



Kurucz (1995) 

Bessell et al. (1989,1991) 

Fluks et al. (1994) 

Tefi 

3500 ~ 50,000 K 

2500 ~ 3800 K 

2500 ~ 3800 K 

log g 

0.0 ~ 5.0 

-1.0 ~ +1.0 

red giant sequence 

[M/H] 

-5.0 ~ +1.0 

-1.0 ~ +0.5 

solar 

AA(nm) 

9.1 ~ 160,000 

491 ~ 4090 

99 ~ 12,500 

n(binA) 

1221 

705 

10,912 



wavelength (nm) 


Fig. 2. B-library giant model spectra {logg = —0.29) for a range in 
metallicity. All the spectra are normalized at A = 817nm. 


spectra - which account for the spectral changes due to varia¬ 
tions in metallicity - have been complemented by the F-library 
spectra at ultraviolet wavelengths. 

This step was achieved with the following procedure: 

— All the F- and B-spectra have first been resampled at the 

(same) uniform grid of wavelengths given by the K-library 
spectra. 

— F-spectra were then recomputed for the effective temper¬ 

atures associated with the B-spectra. This was done by 
interpolation of the F-library sequence of 11 spectra repre¬ 
senting M-giants of types MO through MIO, and using the 
spectral type - Teu scale defined by Fluks et at. 

— Finally, in order to avoid the spurious spikes at AASlOnm 

present in the synthetic B-spectra (Worthey 1994), each B- 
spectrum was combined with the blue part (9.9 ~ 600nm) 
of the corresponding F-spectrum, i.e., having the same Tefi 
and having been rescaled to the B-spectrum flux level at 
AAOOOnm. 

The hybrid spectra created in this manner will hereafter 
be called ‘B-fF-spectra’. 

Of course in this way, completeness in wavelength coverage 
can only be established for solar-abundance models. However, 


extensions of the B-spectra to optical wavelengths down to 
the atmospheric limit at A = 320nm are being worked out 
now and will supersede the corresponding preliminary hybrid 
B-|-F-spectra (Buser et al. 1997). 

3. Comparison with empirical calibrations 

In order to assess the reliability of the synthetic spectra, we 
now compare them to empirical temperature-color calibra¬ 
tions. Depending on the availability and the quality of cali¬ 
bration data, two basic calibration sequences will be used for 
the cooler giants and for the hotter main sequence stars, re¬ 
spectively. 

3.1. Empirical color-temperature relations 

3.1.1. Red giants and supergiants 

Ridgway et al. (1980) derived an empirical temperature-(V-K) 
relation for cool giant stars. This relation is based on stellar di¬ 
ameter and flux measurements, and therefore on the definition 
of the effective temperature: 

( 1 ) 

where fboi is the apparent bolometric flux, a is Boltzmann’s 
constant, and (p is the angular diameter. Hence fboi is almost 
entirely empirical. 

Over the range 5000K ~ 3250K, the Ridgway et al. calibra¬ 
tion was adopted as the effective temperature scale for the V-K 
colors. We derived the color-temperature relations for V-I, J- 
H, H-K, J-K and K-L using the infrared two-color calibrations 
given by Bessell & Brett (1988) (hereafter referred to as BB88). 
For the U-B and B-V colors, we used the color-color relations 
established by Johnson (1966), and Bessell’s (1979) calibration 
was adopted for the (R-I)-Teff relation. 

Because existing calibrations do not go below ~ 3200K, we 
have used both observations and theoretical results published 
by Fluks et al. (1994) in order to construct semi-empirical 
Tefi — color calibrations down to the range 3250K ^ 2500K. 
Synthetic V-K colors computed from their sequence of pho- 
tospheric model spectra provide a very good match to the 
calibration by Ridgway et al., which could thus be extended 
to the range (3767K ~ 2500K) by adapting the theoretical 
Tefr-(V-K) relation from Fluks et alB We then used Fluks 

^The F- models were preferred to the B- models in establish¬ 
ing the semi-empirical Teff-(V-K) relation since they include 
a more accurate calculation of the opacity (the Opacity Sam¬ 
pling technique was used instead of the simple Straight Mean 
method which causes problems when lines saturate). These 
models also incorporate atomic data, as well as new opaci- 
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et al. ’s compilation of (UB)j(VRI)c(JHKLM)ESO observations 
of a large sample of bright M-giant stars in the solar neigh¬ 
bourhood to establish the mean intrinsic colors and standard 
deviations from their estimates of interstellar extinction within 
each photometric band. Finally, with these results, we derived 
mean intrinsic color-color relations - (V-I)-(V-K), (U-B)-(V- 
I), (R-I)-(V-K) and (B-V)-(R-I) -, which allowed us to trans¬ 
late to all these other colors the basic Teff-(V-K) relation 
adopted above for red giants within the temperature range 
3250K ~ 2500K. 

Examples of the adopted fits for the (V-I)-(V-K) and the 
(U-B)-(V-I) sequences are shown in Figs. 3 and 4. 
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2 4 6 8 10 

V-K 

Fig. 3. Adopted (V—I)-(V-K) two-color calibration sequence for 
cool giants. Symbols represent mean values derived from obser¬ 
vations given by Finks et al. (1994). The solid line is a linear 
least-squares fit to the data for V-K > 5. 

For the infrared colors, the photometric data given by Fluks 
et al. are defined on the {JHKL)eso filter system, which is 
different from the filter system defined by BBSS. Using the 
color equations relating the two systems and derived by BBSS, 
transformed JHKL colors from the Fluks et al. data were com¬ 
puted. The resulting (V-K)-(V-J), (V-K)-(V-H), (V-L)-(V- 
J) and (V-H)-(V-J) sequences are well approximated by lin¬ 
ear extrapolations of the two-color relations given by BBSS. 
Furthermore, the model colors derived from the Fluks et al. 
synthetic spectra and the filter responses defined by BBSS also 
agree very well with these extrapolated relations (Fig. 5). We 
therefore chose this method to derive the J-H, H-K, J-K and 
K-L colors over the range 3250K ~ 2500K. However, the un¬ 
certainty implied by the extrapolation to the reddest giants is 
of the order of 0.1 ~ 0.2 mag, indicating that for the coolest 
temperatures near 2500 K the resulting empirical calibration 

ties for the molecules {e.g. H 2 O, TiO, VO), hence providing 
significant a improvement in the V-K synthetic colors at solar 
metallicity (Plez et al. 1992). 



V-I 

Fig. 4. Adopted (U-B)-(V— I) two-color calibration sequence for cool 
giants. Symbols represent the same as in Fig. 3. The solid line is a 
quadratic fit to the data for V-I > 2. 

of the J-H, H-K, J-K, and K-L colors should be improved by 
future observations. 

Table 2 presents the final adopted temperature-color cal¬ 
ibrations for red giants, which was supplemented by a log g 
sequence related to the effective temperature via an evolution¬ 
ary track (IM 0 ) calculated by Schaller et al. (1992). 

3.1.2. Main sequence stars 

To construct empirical Teft-color sequences from 12000 K to 
3600 K for the main sequence stars, we used different cali¬ 
brations: Schmidt-Kaler (1982) was chosen to relate Teft to 
U-B, B-V or R-I, and the two-color relations established by 
FitzGerald (1970), Bessell (1979), and BB88 were then used 
to derive the temperature scales for the remaining colors. This 
procedure should provide color-temperature calibrations with 
uncertainties of < 0.05 mag in color or < 200 K in temperature 
(Buser & Kurucz 1992). 

3.2. Comparison of theoretical and empirical color-temperature 
relations 

3.2.1. Red giants and supergiants 

In order to compare the models to the above color-temperature 
relations for red giants, model spectra were first interpolated 
in the theoretical libraries for appropriate values of surface 
gravity given by the log g-Tefr relation defined by the IMq 
evolutionary track from Schaller et al. (1992). Synthetic colors 
computed from these model spectra are then directly compared 
to the empirical color-temperature relations, as illustrated in 
Fig. 6. 

It is evident that the color differences between equivalent 
models from the K- and the B-|-F-libraries can be as high as 
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Table 2. Adopted semi-empirical Teff-color calibrations and log g values for cool giants. (The notes refer to the method 
employed to derive the semi-empirical colors at a given temperature.) 


Teff 

U-B 

B-V 

V-I 

V-K 

R-I 

J-H 

H-K 

J-K 

K-L 

log g 

4593 

1.017'" 

1.094® 

1.080* 

2.500* 

0.487® 

0.580* 

0.100* 

0.680* 

0.080* 

2.85 

4436 

1.187® 

1.173® 

1.170* 

2.700* 

0.530® 

0.630* 

0.115* 

0.740* 

0.090* 

2.50 

4245 

1.399® 

1.281® 

1.360* 

3.000* 

0.602® 

0.680* 

0.140* 

0.820* 

0.100* 

2.12 

4095 

1.566® 

1.364® 

1.479* 

3.260* 

0.672® 

0.730* 

0.150* 

0.880* 

0.110* 

1.82 

3954 

1.714® 

1.443® 

1.634* 

3.600* 

0.773® 

0.790* 

0.165* 

0.950* 

0.120* 

1.55 

3870 

1.784® 

1.489® 

1.768* 

3.850* 

0.859® 

0.830* 

0.190* 

1.010* 

0.120* 

1.39 

3801 

1.815® 

1.524® 

1.899* 

4.050* 

0.948® 

0.850* 

0.205* 

1.050* 

0.130* 

1.25 

3730 

1.812® 

1.552® 

2.053* 

4.300* 

1.058® 

0.870* 

0.215* 

1.080* 

0.150* 

1.15 

3640 

1.750® 

1.577® 

2.269* 

4.640* 

1.228® 

0.900* 

0.235* 

1.130* 

0.170* 

0.98 

3560 

1.651® 

1.590® 

2.472* 

5.100* 

1.568* 

0.930* 

0.245* 

1.170* 

0.180* 

0.83 

3420 

1.412® 

1.589® 

2.828* 

5.960* 

1.899* 

0.950* 

0.285* 

1.230* 

0.200* 

0.56 

3250 

1.019'‘ 

1.527* 

3.309* 

6.840* 

2.170* 

0.960* 

0.300* 

1.260* 

0.256* 

0.21 

3126 

0.645^ 

1.499* 

3.709* 

7.830® 

2.391* 

0.950^* 

0.370® 

1.320® 

0.310® 

-0.05 

2890 

0.096'‘ 

1.512* 

4.234* 

8.760® 

2.519* 

0.920® 

0.400® 

1.320® 

0.420® 

-0.57 

2667 

-0.146'* 

1.507* 

4.439* 

9.310® 

2.558* 

0.900® 

0.410® 

1.310® 

0.510® 

-1.09 

2500 

-0.328* 

1.510* 

4.593* 

9.560® 

2.567* 

0.880® 

0.420® 

1.300® 

0.550® 

-1.52 


^ Bessell & Brett (1988) two-color relation with Ridgway et al. (1980) to relate Teff to V-K. 
^ Extrapolation of two-color relations from Bessell & Brett (1988). 

^ Synthetic color indices from Finks et al. (1994) models. 

^ From mean two-color relations derived from the Fluks et al. (1994) observed data. 

^ Empirical calibration from Johnson (1966). 

® Empirical calibration from Bessell (1979). 



V-J 

Fig. 5. Adopted extrapolation of the (V—K)-(V-J) relation for cool 
giants (dashed line). Model colors (open squares) and observed data 
(crosses) from Fluks et al. (1994) are compared to the empirical 
sequence adopted. (See text for explanations.) 


0-4 mag, while those between the theoretical library spectra 
and the empirical calibrations may be even larger, up to 1 
mag. 

Such differences - both between the original libraries and 
between these and the empirical calibrations - make it clear 
that direct use of these original theoretical data in population 
and evolutionary synthesis is bound to generate a great deal of 
confusion in the interpretation of results. In particular, appli¬ 
cations to the integrated light of galaxies at faint magnitude 
levels, where effects of cosmological redshift may come into 
play as well, will provide rather limited physical insight unless 
the basic building blocks of the evolutionary synthesis - i.e., 
the stellar spectra - are systematically consistent with the best 
available observational evidence. Thus, our work is driven by 
the systematic consistency of theoretical stellar colors and em¬ 
pirical calibrations as a minimum requirement for the (future) 
standard library. As a viable operational step in this direction, 
a suitable correction procedure for the theoretical spectra will 
be developed in the following section. 

3.2.2. Main sequence stars 

The same procedure as for the giants was applied for the main 
sequence stars, except that a zero-age main sequence isochrone 
(ZAMS) compiled by Bruzual (1995) was used in the appro¬ 
priate interpolation for the surface gravity log g. Again, syn¬ 
thetic photometry results obtained from the theoretical library 
(Kurucz) are compared to the empirical color-temperature re¬ 
lations in Fig. 7. 

Note that the differences between the theoretical and the 
empirical colors are significantly smaller than those for the gi¬ 
ants given in Fig. 6: for the hotter temperatures, they do not 
exceed 0.1 mag., while at cooler temperatures (~ 3500 K) dif- 
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4000 3000 2000 



4000 3000 2000 


4000 3000 2000 



Fig. 6. Empirical color—effective temperature calibrations for solar—metallicity red giant stars (solid lines, according to Table 2) compared 
to the corresponding theoretical relations calculated from original synthetic library spectra (symbols, according to key in insert). Note that 
different scales have been used for the different colors. 


ferences of up to 0.3 mag in V-I between models and obser¬ 
vations again indicate that the coolest K-library spectra still 
carry large uncertainties and should, therefore, be used with 
caution (e.g., Buser & Kurucz 1992). Thus, application of the 
same correction procedure as for the giant models appears war¬ 
ranted for the dwarf models as well. 

Also note that the Kurucz spectra only go down to 3500K. 
We are thus missing the low-luminosity, low-temperature main 
sequence M stars in the present library. However, we anticipate 
here that in a corollary paper (Lejeune et al. 1997, hereafter 
Paper II) the necessary extension is being provided from a 
similar treatment of the comprehensive grid of M-star model 
spectra published by Allard & Hauschildt (1995). 

4. Calibration algorithm for theoretical library spectra 

We now establish a correction procedure for the library spectra 
which preserves their detailed features but modifies their con- 
tinua in such a way that the synthetic colors from the corrected 
library spectra conform to the empirical color-temperature re¬ 
lations. Since the empirical color-temperature relations do not 
provide direct access to the stellar continua, pseudo-continua 
are instead being calculated for each temperature from both 


the empirical (Table 2) and the theoretical (model-generated) 
colors. The ratio between the two pseudo-continua then pro¬ 
vides the desired correction function for the given Tefi. 

4.1. Pseudo-continuum definition 

For a given stellar flux spectrum of effective temperature Te//, 
we define the pseudo-continua pcx{Teff) as black bodies of 
color temperature T'c(A) varying with wavelength: 

pcx{nff) = ainff)-Bxinix)), (2) 

where ce{Teff) is a scale factor and B\{T) is the black body 
function, both of which need to be determined by least-squares 
fits of pc\{Tsff) to the broad-band fluxes of the given flux 
spectrum. Thus, 

j Bx{Tmear,)Si{X)dXc^ fi,i = l,...,9, (3) 

where Si is the normalized transmission function of the 
passband i and fi is the broad-band flux measured through 
this passband. Because colors are relative measurements, we 
normalize by arbitrarily setting the absolute flux in the I-band 
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l.SxlO'* 8000 4000 



1.2x10^ 8000 4000 



1.2x10^ 8000 4000 


Fig. 7. Empirical color—effective temperature calibrations for solar—metallicity dwarf stars (solid lines, see text for sources) compared to 
the corresponding theoretical relations calculated from original synthetic K—library spectra (symbols). Note that different scales have been 
used for the different colors. 


to be equal to 100: fi = 100. The black-body fit in eqn. (3) is 
then obtained iteratively by a conjugate gradients method. 

Fig. 8 illustrates a typical result. Note that, because effects 
of blanketing are ignored by this fitting procedure, the mean 
temperature, Tmean, associated with the best-fitting black- 
body curve, is systematically lower than the effective temper¬ 
ature of the actual flux spectrum. 

a{Teff) having thus been determined, the color tempera¬ 
tures, Tc(Ai), can be derived in a straightforward manner at 
the mean wavelengths □ Ai of the passbands i via the equations: 

a{Teff)-J Bx{n{\i))Si{X)dX = (4) 


®The mean wavelength Ai of a filter of transmission function 
Si is defined in the following way: 

_ f XSi{X)dX 
^ fSi(x)dx ■ 


Interpolation between the Ai by a spline function finally 
provides the continuous (and smooth) color temperatures Tc{X) 
(Fig. 10) required to calculate the pseudo-continua defined by 
equation (2). 


4-2. Correction procedure 

The correction procedure is defined by the following sequence, 
and illustrated (steps 1 to 4) in Fig. 10. 

1. At effective temperature Te//, the empirical pseudo¬ 
continuum (Te//) is computed from the colors of the 
empirical temperature-color relations given in Table 2. 

2. At the same effective temperature T),//, the synthetic 
pseudo-continuum is computed from the synthetic colors 
obtained for the original theoretical solar-abundance spec¬ 
tra given in the K- and/or B-|-F-libraries: 

pcl^"iTeff) =pcl^"{Teff,logg € sequence, \M/H] = 0),(5) 


where log g is defined in the same way as in Sect. 3 by a 
IMq evolutionary track calculated by Schaller et al. (1992). 


















































































5. Finally, the normalization of flnxes in the I-band, ft = 100, 
adopted initially for calculating the pseudo-continua must 
be cancelled now in order to restore the effective temper¬ 
ature scale. Thus, each corrected spectrum of the library, 
fx,corr is Scaled by a constant factor, ^{T,logg,x), to give 
the final corrected spectrum fx^corr- 

fx,corr{T, logg, x) = ^{T, logg, x) ■ fx,corr{T, logg, x), (8) 


where 

aT,logg,x) = J f^(T, logg,x) ■ <^>xiT) ■ dX 


( 9 ) 


assures that the emergent integral flux of the final cor¬ 
rected spectrum conforms to the definition of the effective 
temperature (eqn. (1)). 


The final corrected spectra are thus in a format which al¬ 
lows immediate applications in population evolutionary syn¬ 
thesis. For a stellar model of given mass, metallicity, and age, 
the radius R is determined by calculations of its evolutionary 
track in the theoretical HR diagram, and the total emergent 
flux at each wavelength can hence be obtained from the present 
library spectra via 


Fx{T, logg, x) = 47rR^ • fx,corr(T, logg, x)- (10) 


Fig. 8. Black-body fit of temperature Tmean to a synthetic spec¬ 
trum for Teff = 3126K, illustrating the least-squares solution of 
equations (3). The crosses represent the integrated (heterochro- 
matic) fluxes of the synthetic spectrum measured in the different 
(broad) bandpasses. 


3. The correction function is calculated as the ratio of the 
empirical pseudo-continuum and the synthetic pseudo¬ 
continuum at the effective temperature Tez/u 


^x{T) 


PcT^jT.ff) 
PcTiT.ff) ■ 


( 6 ) 


4. Corrected spectra, carry then calculated from the 
original library spectra, f\-. 


fx.aorr(T,logg,x) = fx(T,logg,x) ■ ^xiT). (7) 


Note that the correction defined in this way becomes an ad¬ 
ditive constant on a logarithmic, or magnitude scale. There¬ 
fore, at each effective temperature the original monochro¬ 
matic magnitude differences between model library spectra 
having different metallicities [M/H] and/or different surface 
gravities log g are conserved after correction. As we shall 
see below, this will be true to good approximation even 
for the heterochromatic broad-band magnitudes and col¬ 
ors because, as shown in Fig. 12, the wavelength-dependent 
correction functions do not, in general, exhibit dramatically 
changing amplitudes within the passbands. 

Fig. 11 shows how the resulting correction functions change 
with decreasing effective temperature. Finally, Figs. 12 
and 13 display the corresponding effective temperature se¬ 
quence of original and corrected spectra at different metal¬ 
licities for the full (Fig. 12) and the visible (Fig. 13) wave¬ 
length ranges, respectively. 


"^Note that in some equations that follow we use a more 
compact notation to designate fundamental stellar parame¬ 
ters which is related to the usual notation by the equivalence 
{T,logg,x) = {Tcff,logg,[M/H]). 


5. Results: the corrected library spectra 

We now discuss the properties of the new library spectra which 
result from the correction algorithm developed above, and 
which are most important in the context of population and 
evolutionary synthesis. 

5.1. Tes-color relations 

Fig. 14 illustrates the Tcff-color relations obtained after cor¬ 
rection of the giant sequence spectra from the K- and B-l-F- 
libraries. Comparison with the corresponding Fig. 6 for the 
uncorrected spectra shows that the original differences which 
existed both between overlapping spectra of the two libraries 
and between the synthetic and empirical relations have indeed 
almost entirely been eliminated. While remaining differences 
between libraries are negligible, those between theoretical and 
empirical relations are below 0.1 mag. 

Fig. 15 illustrates similar results for the main sequence. 
Again, comparison with the corresponding Fig. 7 before correc¬ 
tion shows that the present calibration algorithm provides the¬ 
oretical color-temperature relations which are in almost per¬ 
fect agreement with the empirical data. 

Thus at this point, we can say that for solar abundances, 
the new library provides purely synthetic giant and dwarf 
star spectra that in general fit empirical color-temperature 
calibrations to within better than 0.1 mag over significant 
ranges of wavelengths and temperatures, and even to within 
a few hundredths of a magnitude for the hotter temperatures, 
Tcff > 4000A. 

5.2. Bolometric corrections 

Bolometric corrections, BCv, are indispensable for the direct 
conversion of the theoretical HR diagram, MsoiiTcf f), into the 
observational color-absolute magnitude diagram, Mv{B — V): 


BCv = Mi,oi — Mv + constant. 


( 11 ) 
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Fig. 10. Correction procedure. All fluxes are normalized to fj = 100. 


where the bolometric magnitude, 

poo 

^^bol ~ 2,.3log I fx,corrdX 

Jo ( 12 ) 

= -2.5log{aT^„/7T), 

provides the direct link to the effective temperature (scale). 
Of course, bolometric corrections applying to any other (arbi¬ 
trary) passbands are then consistently calculated from BCi = 
BCv + {Mv — Mi), where the color Mv — Mi is synthesized 
from the corrected library spectra. 

Fig. 16 provides a representative plot of bolometric cor¬ 
rections, BCv, for solar-abundance dwarf model spectra. The 
arbitrary constant in Eq. 11 has been defined in order to 
fix to zero the smallest bolometric correction (Buser & Ku- 
rucz 1978) found for the non-corrected models, which gives 


BCq — 567(5577^4.44^0.0) = —0.190. Comparison with the em¬ 
pirical calibration given by Schmidt-Kaler (1982) demonstrates 
that the present correction algorithm is reliable in this respect, 
too: predictions everywhere agree with the empirical data to 
within ~ 0.05mag - which is excellent. Similar tests for the gi¬ 
ant models also indicate that the correction procedure provides 
theoretical bolometric corrections in better agreement with the 
observations. These results will be discussed in a subsequent 
paper based on a more systematic application to multicolor 
data for cluster and field stars (Lejeune et at 1997). 

5.3. Grid of differential colors 

Since comprehensive empirical calibration data have only been 
available for the full temperature sequences of solar-abundance 
giant and dwarf stars, direct calibration of the present library 
spectra using these data has, by necessity, also been limited to 
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Fig. 12. Normalized corrected (solid lines) and original (dashed lines) library spectra for ranges in effective temperature and metallicity 
and covering wavelengths from the photometric U- through K-passbands. Top panels: K—library dwarf models; middle panels: K-library 
giant models; bottom panels: B+F-library giant models. 


solar-abundance models. However, because one of the principal 
purposes of the present work has been to make available the¬ 
oretical flux spectra covering a wide range in metallicities, it 
is important that the present calibration for solar-abundance 
models be propagated consistently into the remaining library 
spectra for parameter values ranging outside those represented 
by the calibration sequences. We thus have designed our cor¬ 
rection algorithm in such a way as to preserve, at eaeh temper¬ 
ature, the monochromatic flux ratios between the original spec¬ 


tra for different metallieities [M/HJ and/or surface gravities 
log g. Justification of this procedure comes from the fact that, 
if used differentially, most modern grids of model-atmosphere 
spectra come close to reproducing observed stellar properties 
with relatively high systematic accuracy over wide ranges in 
physical parameters (e.g., Buser & Kurucz 1992, Lejeune & 
Buser 1996). 

In order to check the extent to which preservation of 
monochromatic fluxes propagates into the broad-band colors. 
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Fig. 13. Same as Fig. 12, but for the visible—near infrared wavelength ranges only. 


we have calculated the differential colors due to metallicity dif¬ 
ferences between models of the same effective temperature and 
surface gravity: 

^(Cj,[M/ir]) = = 1, •••, 8. (13) 

We can then calculate the residual color differences between 
the corrected and the original grids: 

5(A(c,,[m/«])) = A(c“[m/«]) - (14) 


Results are presented in Figs. 17 and 18 for the coolest 
K-library models (SSOOTf < Te// < SOOOR') and for the B+F- 
library models for M giants (25007^ <Teff < 3750K), respec¬ 
tively. Residuals are plotted as a function of the model number, 
which increases with both surface gravity and effective temper¬ 
ature, as given in the calibration sequences. The different lines 
represent different metallicities, —3.0 < [M/H] < -1-1.0, as ex¬ 
plained in the captions. 

The most important conclusion is that, in general, the cor¬ 
rection algorithm does not alter the original differential grid 
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Fig. 14. Empirical color—effective temperature calibrations for solar—metallicity red giant stars (solid lines, according to Table 2) compared 
to the corresponding theoretical relations calculated from corrected synthetic library spectra (symbols, according to key in insert). Compare 
with Fig. 6. 


properties significantly for most colors and most temperatures 
- in fact, the residuals are smaller than only a few hundredths 
of a magnitude. Typically, the largest residuals are found for 
the coolest temperatures {T^ff < SSOOif) and the shortest- 
wavelength colors, UBVRI, where the correction functions of 
Fig. 11 show the largest variations not only between the differ¬ 
ent passbands, but also within the individual passbands. This 
changes their effective wavelengths and, hence, the baselines 
defining the color scales (cf. Buser 1978). Since this effect tends 
to grow with the width of the passband, it is mainly the co¬ 
incidence of large changes in both amplitudes and slopes of 
the correction functions with the wide-winged R-band which 
causes residuals for the R-I colors to be relatively large in Fig. 
18. 

Calculations of color effects induced by surface gravity 
changes lead to similar results. This corroborates our con¬ 
clusion that the present correction algorithm indeed provides 
a new model spectra library which essentially incorporates, 
to within useful accuracy for the purpose, the currently best 
knowledge of fundamental stellar properties: a full-range color- 
calibration in terms of empirical effective temperatures at solar 
abundances (where comprehensive calibration data exist) and 
a systematic grid of differential colors predicted by the original 


theoretical model-atmosphere calculations for the full ranges 
of metallicities and surface gravities (where empirical data are 
still too scarce to allow comprehensive grid calibration). 

Of course, we are aware that the present correction algo¬ 
rithm becomes increasingly inadequate with the complexity of 
the stellar spectra growing with decreasing temperature and/or 
increasing surface gravity and metallicity. For example, be¬ 
cause under these conditions the highly nonlinear effects of 
blanketing due to line saturation and crowding and broad¬ 
band molecular absorption tend to dominate the behavior of 
stellar colors, particularly at shorter (i.e., visible) wavelengths, 
even the corresponding differential colors cannot either be re¬ 
covered in a physically consistent manner by a simple linear 
model such as the present. However, the limits of this approach 
will be further explored in Paper II, where the calibration of 
theoretical spectra for M-dwarfs will be attempted by intro¬ 
ducing the conservation of original differential colors of grid 
spectra as a constraint imposed to the correction algorithm. 

6. Organization of the library 

The corrected spectra have finally been composed into the uni¬ 
fied library shown in Fig. 19. In this library, model spectra are 
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Fig. 15. Empirical color-effective temperature calibrations for solar—metallicity dwarf stars (solid lines, see text for sources) compared to 
the corresponding theoretical relations calculated from corrected synthetic K-library spectra (symbols). Compare with Fig. 7. 


given for a parameter grid which is uniformly sampled in Tsff, 
log g, and [M/H], each spectrum being available for the same 
wavelength grid, (9lA < A < 1, 600, OOOA), with a medium res¬ 
olution of 10-20 A. 

The main body consists of the K-library, which provides 
the most extensive coverage of parameter space (cf. Table 1) 
down to Teff = 3500K, and which has therefore been fully 
implemented to this limit. The extensions to lower tempera¬ 
tures, 3500it' > Tsff > 2500it', and associated lower surface 
gravities, 0.87 > logg > —1.02, are provided by the M-giant 
spectra from the B+F-library, which, however, covers only the 
limited metallicity range -1-0.5 > [M/H] > —1.0. 

This unified library is available (in electronic form) at the 
CDS, Strasbourg, France, where it can be obtained in either 
of two versions, with the corrected or the original uncorreeted 
spectra from the K- and B+F-libraries. 

7. Discussion and conclusion 

Although astronomers have for a long time agreed that a uni¬ 
form, complete, and realistic stellar library is urgently needed, 
it must be emphasized that this goal has remained too ambi¬ 
tious to be achieved in a single concerted effort through the 
present epoch. We have thus attempted to proceed in well- 


defined steps, with priorities set according to the availability 
of basic data and following the most obvious scientific questions 
that would likely become more tractable, or even answerable. 
Therefore, we briefly review the present achievement to clar¬ 
ify its status in the ongoing process toward a future standard 
library of theoretical stellar spectra for photometric evolution¬ 
ary synthesis. 

1. Completeness. The unification of the massive K-library 
spectra with those for high- luminosity M-star models (the 
B+F-library) is most important, because even as a small 
minority of the number population of a given stellar system, 
the late-type giants and supergiants may provide a large 
fraction of this system’s integrated light at visible and in¬ 
frared wavelengths. This fact was recognized early on (e.g., 
Baum 1959), and eventually also co-motivated the effort 
leading to the existence of the B-library used in this work 
(Bessell et al. 1988). 

While the K-, B-, and F-libraries have been used to rem¬ 
edy incompleteness (in either wavelength or parameter cov¬ 
erage, or both) of available observed stellar libraries be¬ 
fore (e.g., Worthey 1992, 1994; Fioc & Rocca-Volmerange 
1996), our first goal here has been to join them as a purely 
theoretical library, providing the main advantages of phys- 
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ical homogeneity and definition in terms of fundamental 
stellar parameters - which allows direct use with stellar 
evolutionary calculations. 

But even so, the present library remains incomplete in 
several respects. First, stellar evolution calculations (e.g.. 
Green et al. 1987) predict that high-luminosity stars with 
temperatures near or below 3500Jf may also exist at low 
metallicities, \M/H\ ~ —1.0, and their flux contributions 
at visible-near ultraviolet wavelengths (where metallicity 
produces significant effects) may not quite be negligible in 
the integrated light of old stellar populations. Therefore, 
in order to provide the fuller coverage required for an ade¬ 
quate study of this metallicity-sensitive domain, new calcu¬ 
lations of B-library spectra extending the original data to 
both shorter wavelengths (AA > 320nm) and lower abun¬ 
dances {[M/H\ > —2) (Buser et al. 1997) will replace the 
current hybrid B-|-F spectra and make the next library ver¬ 
sion more homogeneous. 

Secondly, even though the low-temperature, 
low-luminosity M-dwarf stars do not contribute signifi¬ 
cantly to the integrated bolometric flux, they are not negli¬ 
gible in the determination of mass-to-light ratios in stellar 
populations. Thus, a suitable grid of (theoretical) M-dwarf 
spectra calculated by Allard & Hauschildt (1995) is being 
subjected to a similar calibration process (Paper II) and 
will be implemented in the present library as an important 
step toward the intended standard stellar library. 

Finally, the libraries of synthetic spectra for hot O- and 
WR-stars which were recently calculated by Schmutz et 
al. (1992) will allow us to extend the calibration algorithm 
to ultraviolet lUE colors, where such stars radiate most of 
their light. 

2. Realism. In view of its major intended application - pho¬ 
tometric evolutionary synthesis -, the minimum require¬ 
ment that we insist the theoretical library must satisfy, is to 
provide stellar flux spectra having (synthetic) colors which 
are systematically consistent with calibrations derived from 
observations. How else could we hope to learn the physics 
of distant stellar populations from their integrated colors, 
unless the basic building blocks - i.e., the library spec¬ 
tra used in the synthesis calculations - can be taken as 
adequate representations of the better-known fundamen¬ 
tal stellar properties, such as their color-temperature rela¬ 
tions? 

Because the original library spectra do not meet the above 
minimum requirement (Sect. 3), we have developed an al¬ 
gorithm for calibrating existing theoretical spectra against 
empirical color-temperature relations (Sect. 4). Because 
comprehensive empirical data are unavailable for large seg¬ 
ments of the parameter space covered by the theoretical li¬ 
brary, direct calibration can be effected only for the major 
sequences of solar-abundance models (Sect. 5). However, 
we have also shown that the present algorithm provides 
the desired broad-band (or pseudo-continuum) color cali¬ 
bration without destroying the original relative monochro¬ 
matic fluxes between arbitrary model grid spectra and 
solar-abundance calibration sequence spectra of the same 
effective temperature. This conservation of original grid 
properties also propagates with useful systematic accu¬ 
racy even through most differential broad-band colors of 
the corrected library spectra. Thus, to the extent that 


differential broad-band colors of original library spectra 
were previously shown to be consistent with spectroscopic 
or other empirical calibrations of the UBV-, RGU-, and 
Washington ultraviolet-excess-metallicity relations (Buser 
& Fenkart 1990, Buser & Kurucz 1992, Lejeune & Buser 
1996), the corrected library spectra are still consistent with 
the same calibrations. 

3. Library development. At this point, we feel that some 
of the more important intrinsic properties required of the 
future standard library have already been established. Of 
course, many more consistency tests and calibrations will 
now be needed that can, however, only be performed for 
local volumes of the full parameter space covered by the 
new library. For example, we shall use libraries of observed 
flux spectra for individual field and cluster stars to better 
assess - and/or improve - the performance of the present 
library version in the non-solar-abundance and non-visual 
wavelength regimes. Eventually, we also expect significant 
guidance toward a more systematically realistic version of 
the library from actual evolutionary synthesis calculations 
of the integrated spectra and colors of globular clusters 
(Lejeune 1997). 

Last, but not least, we would like to emphasize that, while 
we here present the results taylored according to the gen¬ 
eral needs in the field of photometric evolutionary synthe¬ 
sis, the library construction algorithm has been designed 
such as to allow flexible adaptation to alternative calibra¬ 
tion data as well. As we shall ourselves peruse this flexi¬ 
bility to accommodate both feed-back and new data, the 
reader, too, is invited to define his or her own preferred 
calibration constraints and have the algorithm adapted to 
perform accordingly. 

Acknowledgements. We are grateful to Michael Scholz and 
Gustavo Bruzual for providing vital input and critical discus¬ 
sions. Christophe Pichon is also aknowledged for his precious 
help with the final implementation of some of the figures. We 
wish to thank warmly the referee for his helpful comments and 
suggestions. This work was supported by the Swiss National 
Science Foundation. 


References 

Allard F., Hauschildt P.H., 1995, ApJ, 445, 433 
Baum W. A., 1959, PASP, 71, 106 

Bell R.A., Eriksson K., Gustafsson B., Nordlund A., 1976, 
A&AS 23, 37 

Bessell M.S., 1979, PASP, 91, 589 
Bessell M.S., Brett J.M., 1988, PASP, 100, 1134 
Bessell M.S., Brett J.M., Scholz M., Wood P.R., 1989, A&AS, 
77, 1 

Bessell M.S., Brett J.M., Scholz M., Wood P.R., 1991, A&AS, 
89, 335 

Bruzual G.A., 1995, private communication 
Bruzual G.A., Chariot S., 1993, ApJ, 405, 538 
Bruzual G.A. et al., 1997 (in preparation) 

Buser R., 1978, A&A, 62, 411 

Buser R., Eenkart, R.P. 1990, A&A 239, 243 

Buser R., Kurucz R., 1978, A&A, 70, 555 

Buser R., Kurucz R., 1992, A&A, 264, 557 

Buser R., Scholz M., Lejeune Th. 1997, A&A (in preparation) 



15 


Eriksson K., Bell R.A., Gustafsson B., Nordlund A., 1979, 
Trans. lAU Vol. 17A, Part 2, Reidel, Dordrecht, p. 200 
Fioc M., Rocca-Volmerange B., 1996, A&A, submitted 
FitzGerald M.P., 1970, A&A, 4, 234 
Fluks M.A. et ai, 1994, A&AS, 105, 311 
Green E.M., Demarque P., King G.R., 1987, The Revised Yale 
Isochrones and Luminosity Funetions, Yale Univ. Obs., 
New Haven, USA 

Gustafsson B., Bell R.A., Eriksson K., Nordlund A., 1975, 
A&A 42, 407 

Johnson H.L., 1966, ARA&A 4, 193 
Kurucz R.L., 1979a, ApJS, 40, 1 

Kurucz R.L., 1979b, in: Problems of Calibration of Multicolor 
Photometric Systems, ed. A.G. Davis Philip, Dudley Obs. 
Rep. 14, p. 363 

Kurucz R.L., 1992, in: The Stellar Populations of Galaxies, eds. 
B. Barbuy, A. Renzini, lAU Symp. 149, Dordrecht, Kluwer, 
p. 225 

Kurucz R., 1995, GD-ROM, private communication 
Kurucz R.L., Peytremann E., 1975, Smithsonian Astrophys. 

Obs. Special Rep. 362 
Leitherer G. et ai, 1996, PASP 108, 996 
Lejeune Th., 1997, Ph.D. thesis, in preparation 
Lejeune Th., Buser R., 1996, Baltic Astronomy, 5, 399 
Lejeune Th., Buser R., Guisinier F. 1997, A&A (in preparation. 
Paper II) 

Lejeune Th., Lastennet E., Valls-Gabaud D, Buser R. 1997, 
A&A (in preparation) 

Plez, B., Brett, J., Nordlund, A., 1992, A&A 256, 551 
Ridgway S.T., Joyce R.R., White N.M., Wing R.F., 1980, ApJ, 
235, 126 

Schaller G., Schaerer D., Meynet G., Maeder A., 1992, A&AS, 
96, 269 

Schmidt-Kaler, Th., 1982, in: Landolt-Bor stein, Neue Serie, 
Gruppe VI, Bd. 2b, eds. K Schaifers, H.H. Voigt, Springer, 
Berlin Heidelberg New York, p. 14 
Schmutz W., Leitherer G., Gruenwald R., 1992, PASP, 104, 
682 

Worthey G., 1992, PhD thesis, Univ. of Galifornia, Santa Cruz 
Worthey G., 1994, ApJS, 95, 107 



wavelength (nm) 

Fig. 9. Color temperatures for synthetic spectra covering a range of 
effective temperatures, as labelled. Top pair: K-library dwarf mod¬ 
els; middle pair: K—library giant models; bottom pair: B-l-F-library 
giant models. Note that systematically Tc <Teff, as tracked down 
by the black-body fit temperatures, Tmean- 
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Fig. 11. Correction functions for a range of effective temperatures. 
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Fig. 16. Bolometric correction, BCy, as a function of effective 
temperature, for corrected K-library solar—metallicity dwarf model 
spectra. The solid line represents the empirical calibration given by 
Schmidt-Kaler (1982). 
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Fig. 17. Residuals between differential colors obtained from orig¬ 
inal and corrected library spectra for models of the same effective 
temperatures and surface gravities but different metallicities. The 
abscissa represents the monotonic increase in temperature (as indi¬ 
cated along the top of the panel) and surface gravity in a similar way 
as Table 2. Residuals are for the coolest models from the K-library 
(3500K ~ 5000K) and for metallicities [M/H] = —3.0 (solid line), 
[M/H] = —1.0 (long-dashed line), and [M/H] = -^1.0 (dotted line). 
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Fig. 19. The (Teff, logg, [M/H]) parameter space covered by the 
present unified library. Solid lines delineate approximate boundaries 
of available models. (B+F)-models for M—giants are confined to the 
metallicity range —1.0 < [M/H] < +0.5 (grey part). 


Fig. 18. Same as Fig. 17, but for the M-giant models from the 
B+K—library (2500iC 3750iC) and for metallicities [M/H] = —1.0 
(solid line) and [M/H] = +0.5 (dotted line). The zero residuals for 
U — B are due to the fact that in this wavelength range, B+K-library 
spectra are available for solar metallicity only. 
































































